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We investigate the symmetry of the superconducting order parameter in the proximity of a 
phase-separation or of an incommensurate charge-density-wave instability. The attractive effec- 
tive interaction at small or intermediate transferred momenta is singular near the instability. This 
strongly g-dependent interaction, together with a residual local repulsion between the quasiparticles 
and an enhanced density of states for band structures appropriate for the high temperature super- 
conducting oxides, strongly favors the formation of d-wave superconductivity. The relative stability 
with respect to superconductivity in the s-wave channel is discussed in detail, finding this latter 
hardly realized in the above conditions. The superconducting temperature is mostly determined 
by the closeness to the quantum critical point associated to the charge instability and displays a 
stronger dependence on doping with respect to the simple proximity to a Van Hove singularity. 
The relevance of this scenario and the generic agreement of the resulting phase diagram with the 
properties displayed by high temperature superconducting oxides is discussed. 
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I. INTRODUCTION 



An increasing complexity in the structure of the or- 
der parameter of the high Tc superconductors is coming 
out from the experiments. Strong anisotropic behaviour 
in BiSCCO materials is evident from the photoemission 
experiments. A very small value of the gap is measured 
along the TY direction (| kx \ — \ ky \) compatible with a 



line of nodes according to the d-wave pairingd'EI. Evidence 
for d-wBMe pairing comes also from penetration depth 
measurements and from several Josephson-coupling ex- 
periments mostly in YBCOa. Phase sensitive experi- 
ments, however, in, some cases also provide evidence for 
the s-wave pairingp. The actual situation seems therefore 
to be more complex than the one associated to a specific 
pairing of a given symmetry. Even though the present 
experimental situation is still uncertain it is worth con- 
sidering the implications as well as the various possible 
origins of superconducting order parameters of different 
symmetries. 

In this paper we analyse the symmetry (d-wave vs s- 
wave) and the doping dependence of the critical tem- 
perature Tc for the pairing interactions which arise near 
a charge instability. Specifically, we will show that a 
d-wave pairing comes out in a region of the phase dia- 
gram near a phase-separation (PS) or an incommensu- 
rate charge density wave (CDW) instability. The lat- 
ter situation arises in strongly correlated systems when 
Coulomb forces forbid a thermodynamic PS giving rise 
to a density instability at a specific finite wave vector 
qcDW- In a recent letter it was shown that near these 
charge instabilities a singular scattering among quasipar- 
ticles arises which may be responsible foe the anomalous 
behavior of the normal phase above T^. According to 
our analysis, the same singular scattering may provide 



the strong pairing mechanism needed for high Tc cuprate 
superconductors. The order parameter turns out to be 
strongly anisotropic and, under quite general conditions, 
of d-wave symmetry. This happens in the proximity of 
both PS and CDW. In addition, the analysis of the CDW 
case reveals an interesting interplay between the finite- 
momentum instability and the geometry of the Fermi sur- 
face leading to the appearence of specific structures in 
the superconducting gap. Finally, by considering the PS 
and the CDW instabilities at zero temperature within the 
context of quantum transitions, we will argue on a possi- 
ble mechanism for the existence of both strong variations 
and plateaus in the dependence of the superconducting 
critical temperature on doping. 



II. PROXIMITY TO PS 

The origin of d-wave pairing is usually attributed to the 
relevance of electron-electron (e — e) scattering at large 
wave vector, specifically at q =r(K, tt) due to strong anti- 
ferromagnetic spin fluctuationsflfl. Alternatively, charge 
fluctuations near a PS (in particular driven by excitonic 
effects in a three-band extended Hubbard model) i^*^ 
already been proposed as a source of d-wave pairingSilil. 
Here we shall exploit the general feature that an effectivgj 
interaction with on site repulsion and singular attractionEl 
at small q is generated nea*l«2-|a PS, regardless of the 
forces driving the instabilityOil3. 

An effective e — e interaction with attraction at small 
q and weak on site repulsion was already advocated 
to explain the strong amaatropy in the gap observed 
in ARPES experimental^"!!^. In these analyses the 
anisotropy of the gap (of s-wave symmetry) is driven 
by the anisotropy of the density of states due to mo- 
mentum decoupling induced by the small q scattering. 



1 



which for a given point k in momentum space couples 
only nearby states. Later the analysis was extended by 
considering d-wave pairing which indeed turns out to be 
favored when the local repulsion is sizableEZl. Our follow- 
ing analysis of the symmetry of the order parameter in a 
two-dimensional system naar a PS shares many features 
with the analyjSis of RefO, eventhough the model con- 
sidered in RefO assumes that the small-q attractiaa is 
of pure phononic origin with no connection with PSiia. 

In the proximity of PS the effective interaction has a 
rich dynamical structuije which is relevant for the anoma- 
lous behavior above Tj3. Concerning the superconduct- 
ing properties we limit ourselves to consider the static 
part of this effective interaction which has the formcl 



V 



(1) 



where k is the cutoff of the attractive interaction, V the 
attraction strenght in unity of the inverse square of the 
lattice constant and U is an on site effective repulsiorO; 
we impose to the interaction the periodicity of the lat- 
tice without modifing the behaviour at small q by writing 
ql +qy as 2(2 — {cos{qx) +cos{qy)). Near PS vanishes 
as S ^ S^^ where S is the doping (with respect to half 
filling) and S^^ is the critical doping at which the insta- 
bility for PS occurs. At 6^^ the value V stays finite and 
a singular effective attraction at small q arises. 

To study the symmetry of the superconducting gap we 
solve a simple BCS self-consistent equation for the gap 
parameter A(fc) 



A(fc) 



N ^ 

p 



^ tanh . , ^ 



(2) 



with Veff = Vjff given in Eq.(|]). Later we shall use 
a different effective potential {— V^^^ , see following 
Eq.(^) to analyse the gap in the proximity of a CDW. 

Of course, approaching d^^ the superconducting prop- 
erties obtained within the BCS theory should only be 
considered on a qualitative ground. For instance, wave 
function corrections become relevant. (The same com- 
ment would apply for 5 — > 5^^^ , where 5^^^ is the 
critical doping for the CDW instability). However, it is 
worth noting that vertex corrections could partially com- 
pensate wave function corrections due to small-g iScat ter- 
ing similarly to whaL,happens in one dimensionEiiEil, or 
near one dimensionEj. An explicit calculation has in- 
deed shown that the superconducting critical tempera- 
ture is enhanced by considering the vertex corrections to 
Migdal-Eliashberg theor*-|in a model with small-q cou- 
pling to optical phononsE3. 

In Eq.(|^) the sum over pis done over the first Brillouin 
zone, N is the number of lattice sites and ^1 = ?| + ^| 
with ^j: being the electronic dispersion measured with 
respect to the Fermi energy Ep. We consider a tight- 
binding model with hopping up to the fifth nearest neigh- 
bors 



^(fca;, ky) ^ ^ CiTli(kx-I ky) 



(3) 



where, according to Ref. 
ci = 0.1305eV, f]2{kx^k, 



V) 



-0.5951eV, mikx,ky, 
m{kx,ky) = i(cos2fc: 
Vbikxyky) = ^{cos2kx COS ky 



we choose rji{kx,n,y 

= i(cOs(fca, 

cos fcr COS k. 



cos 2ky) 



ky) = 1 , 

f COs(fcj^)), C2 = 

y, C3 0.1636eV, 
C4 = -0.0519eV, 



COS 2fcy cos /cj;), C5 



— O.lllTeV, riQ{krc,ky) = cos 2A:a; cos 2A:y, and cg = 
0.0510eV. These parameters are appropriate for the band 
structure of the BiSCCO compounds, thus giving an open 
Fermi surface and a van Hove singularity (VHS) slightly 
below the Fermi level (for electrons). The value of Ep = 
— ci ~ — 0.1305eV is fixed to get the proper distance of 
the Fermi surface from the VHS {Ep — Evhs = 35meV 
as determined experimentally), and corresponds to (op- 
timal) doping 5 = 0.17. The fuU bandwidth W is 1.4eV. 

The choice of the parameters entering the effective in- 
teraction Vj'A is to some extent arbitrary since they will 
depend on the specific mechanism driving PS. We choose 
V = Q.21W = 0.3ey and U = Q.2eV. These values are 
of the order of magnitude of the estimates obijainpd for a 
single-band Hubbard-Holstein model near PSE3,EI. 

The BCS equation (^ is solved in two limiting cases: 
i) at zero temperature, to obtain the ground state in 
the superconducting phase; M)in the linearized form, for 
A(fc) going to zero, to evaluate the critical temperature 
Tc of the superconducting transition. We provide a nu- 
merical self-consistent solution of the BCS equation on 
a lattice of = 128 x 128 points (in each quadrant of 
the BrillouiBj-Zone), taking advantage of the Fast Fourier 
TransformnScj. In the case of an undistorted tetragonal 
lattice, both in the proximity to PS and to CDW, the 
solutions in the s-wave and d-wave channels are decou- 
pled and the solution is a superposition of all harmon- 
ics of a given symmetry. In the range of parameters we 
have considered the s-wave solutions can be roughly ap- 
proximated by a linear combination of the first square 
harmonics 



As{k) ~ Ao 4- Ai(cos k^ + cos ky) + A2 cos k^ cos ky 



(4) 



while the d-wave solutions are decoupled in the (i^2_„2- 



like and dxy-^ikc harmonics 



^rf,2_„2 (fc) - A.i{coskx - cosky) 
'^d^yik) ~ Ai sinkx sinky + ... 



(5) 
(6) 



In all cases we find that the dxy-\i^e solution is strongly 
suppressed because it has a line of nodes along the rel- 
evant TM directions where saddle points are present at 
(±7r, 0), (0, ±7r), giviag rise to Van Hove singularities in 
the density of state£3. 

We now proceed to analyse the gap parameter at 
T = OK in the proximity of PS. The solution is given 
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for the same set of parameters at varius k's. In Fig. la 
we show the momentum dependence of the s-wave solu- 
tion on the Fermi surface. We define an angular variable 
(j) measured from the line AIY to detect the points on 
the Fermi surface as seen from the point Y = (tt, tt). The 
s-wave solution has two nodes in each quadrant, sym- 
metrically located with respect to the TV direction (as 
in Ref.Lj), whose positions depend on k. The maximum 
value of the gap is in the FM directions, near the van 
Hove singularities (this is due to almost complete mo- 
mentum decoupling at small q). 
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Fig.l: Angular dependence of the order parameter nor- 
malized by its maximum value for various cutoffs k of the 
attractive interaction close to PS [Eq. (jl|)]. U = 0.2eV 
and V — 0.3eV. The angle is defined in the inset, where 
the Fermi surface branch in the upper right quadrant of 
the Brillouin zone is shown, (a) Angular dependence of the 
s-wave order parameter normalized by its maximum value 
AIj^^ = 237.9, 161.1, 9.4 iC for k = ^, ^, f respectively; 
(b) Angular dependence of the d-wave normalized order pa- 



rameter with Afi^^ = 358.6, 308.0, 85.: 
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respectively. The pure d-wave solution (dot-dashed line) and 
a weaker couplig {V = O.leV) (long-dashed line to be dis- 
cussed later) solution are also shown. 

In Fig. lb we show the same ^-dependence for the 



d^2_y2 solution. This solution has a line of nodes along 
the TY direction (0 = 45°) and the maximum value along 
FM directions. The maximum value of the gap param- 
eter is much larger for the d-wave solution than for the 
s-wave solution and the relative difference is so neat that 
it is easy to predict also a large difference in the conden- 
sation energy. 

The comparison between s-wave and c?-wave is com- 
pleted evaluating the condensation energy per particle 
6F = Fsuper — Fnormai for the BCS grouud state. In 
Fig. 2 we report the results for the two analysed symme- 
tries and different k. The difference between the two con- 
densation energies is considerable, since we find jp- > 3.5 
for all considered k's. 
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Fig.2: 



(solid line) and d-wave (dashed line) conden- 



sation energy per particle SF = Fs, 
of K. U = 0.2eV and V = 0.3eV. 



-F„ 



as a function 



From the above analysis it is clear that the d-wave 
solution is the favorite gap parameter of the supercon- 
ducting phase induced by the interaction (^) character- 
istic of models that present a small-momentum charge 
instability. We would like to point out that this result 
is a consequence of three related but distinct features of 
the system: first there is a strongly g-dependent attrac- 
tion peaked at small momenta. This leads to a momen- 
tum decoupling so that the modulus of the gap parameter 
tends to match the local density of states of the system. 
Then a second feature needed to obtain a gap anisotropy 
is a largely anisotropic density of states, possibly with 
extended van Hove singularities in some regions of the 
k-space. Both these features, however, only point to- 
wards the setting in of an anisotropic order parameter, 
which still could well be of the s type. It is the presence 
of a (nearly) momentum independent sizable repulsion 
which leads to d-wave symmetry. This solution is able to 
avoid the isotropic repulsive interaction U since its av- 
erage over the Brillouin zone is zero while keeping the 
paired electrons in the attractive region of the effective 
potential. Roughly, the d-wave becomes favorable when 
the average repulsion felt by the s-wave paired electrons 
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exceeds the loss in condensation energy due to the van- 
ishing of the order parameter along the nodal regions. 
Among the d waves, the dx2_y2 is preferred because the 
nodes occur where the modulus of the order parameter is 
anyway small (i.e., due to momentum decoupling, where 
the density of states is small). As extensively discussed 
in Appendix A, by reducing the relative value of U with 
respect to V the s-wave solution may be recovered. At 
K — 7r/16 we find that this happens for U* — O.OGeV. We 
believe, however, that the presence of a sizable U > U* 
is an unavoidable feature of strongly correlated systems. 
In turn strong correlation is a prerequisite for the likely 
occurrence of a charge instability and for the related en- 
hancement of the residual attraction among quasiparti- 
cles. This leads to the conclusion that pairing near PS 
in strongly correlated systems is of d-wave symmetry. 



Ix)'^ ~^ ily ~ ly)'^) for (7 q" while mantaining the lattice 
periodicity. 

As in Eq.(|l|), k acts as a cutoff of the effective inter- 
action. It can be identified with the inverse correlation 
length of the CDW and vanishes at 4n^^- 
U and V the same values used beforeo and solve Eq.(||) 
with Veff = Vfjf^ given in Eq.(0). 

All the qualitative features presented above for the case 
of proximity to PS are not changed near the CDW in- 
stability if qcDW is small, specifically if it stays smaller 
than the momentum qp connecting the two branches of 
the Fermi surface around the VHS. For the considered 
set of band parameters we find q*p = 1.2. We report in 
Figs. 3, the values of the gap on the Fermi surface as a 
function of the angle (j) for qcDW = 0-9 both for s and 
d-wave symmetry. 



III. PROXIMITY TO CDW 



The occurence of PS is a theoretical outcome of mod- 
els with short-range forces. In a real charged system 
long-range forces prevent charge segregation on a macro- 
scopic scale leaving open the way to charge segregation 
on a microscopic scale, i.e. to CDWEj. If and when 
CDW are realized is a debated issue. In the early stage 
of the investigation on this problem, a charge-glass be- 
havior was suggested as a result of frustated PSEa. Re- 
cently, CDW have been observed in the related com- 
pound La2-x-yNdySrxC«|Oj23. Superstructures are also 
seen in various cuprateS'E3. It is not yet clear whether 
these experimental evidences can indeed be attribuited 
to the interplay between PS and long range Coulomb 
forces. Here we assume that the system is near a CDW 
instability characterized by an incommensurate wave vec- 
tor qcDW- The value of this qcDW is mainly fixed by 
the balance between charge segregation favored by short 
range forces and the consequent electrostatic cost. To a 
large extent, qcDW is not a Fermi surface property (i.e. 
QCDW is not a nesting vector). In the analysis of the 
Hubbard-Holstein model qcDW is found almost parallel 
to the |(4:7i', 0), (0, ±7r) directions (similarly to the finding 
of Ref.E3 in the analysis of an effective Ising model with 
long range forces). Its value comes out to be of order 
one. 

Near the CDW instability and neglecting the dynam- 
ics relevant to the anomalous properties above Tc, the 
efl^ctive interaction among quasiparticles can be written 



V 



(7) 



where the sum is over the four equivalent vectors of 
the CDW instability q" = {±qcDW,Q),{0,^qcDw) and 
uj^ = 2(2 — cos(g:r - Qx) ^ cos{qy - qy))- This expres- 
sion is used to reproduce the behavior — 1/(k^ + (fe — 
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Fig. 3: Angular dependence of the order parameter nor- 
malized by its maximum value for various cutoffs k of the 
attractive interaction close to a CDW instability [Eq. (§)]. 
U — 0.2eV and V = Q.3eV. (a) s-wave order parameter nor- 
malized by its maximum value A%i^^ = 17.5, 5.9, 1.3 K for 
K = respectively; (b) d-wave normalized order pa- 



rameter, with A^,/„j. 



200.0, 136.4, 21.5 K ior k ^ 



respectively. The pure d-wave solution (dot-dashed line) and 
a weaker couplig {V = 0.2eV and 1^ ~ fg) (long-dashed line) 
solution are also shown. 
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Notice that for all qcDW < Qp the shape of the d-wave 
gap is not very different from a pure dx2-y2 (longdashed 
line in Fig. 3b). The situation would be different if we 
would consider a much smaller value for the attraction 
V (say V ~ O.leV), pushing the system in a weak cou- 
pling regime. In this case there is a tendency of the curve 
to becomCpflat around (j) = tt/A analogously to the find- 
ing of RefJlZI (longdashed curve in Fig. 3b and even more 
evidently, for V ~ O.leV in Fig. lb) 

For qcDW > Qp a. new feature appears in A^. As 
shown in Fig. 4, where we report the gap for various values 
oi qcDW^ by increasing qcDW extended flat regions de- 
velop around the maxima a,t (f) — 0, 7r/2, which eventually 
become local minima for the d-wave solution. This fea- 
ture is very pronounced for small k (k = 7r/16 in Fig. 4). 
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order parameter in the proximity of a CDW instabilities with 
different values of qcDW- U — Q.2eV and V — 0.3eV; (b) 
Same as (a) for d-wave symmetry. 

The reason is the interplay between the VHS, which 
mostly weight 4> = 0,7r/2, and the value of qcDW which 
for qcDW > q*F only connects the two branches of the 
Fermi surface around the VHS at finite values of (j) re- 
sulting in maxima for Afc near these values of (p. Notice 
that these local maxima in the order parameter are much 



more pronounced in the s-wave (Fig. 4a) than in the d- 
wave channel (Fig. 4b) . Moreover, in the s-wave order pa- 
rameter they are also visible when qcDW < q*F provided 
the mass term is large enough [k > 7r/12). In this lat- 
ter case, although the instability wavevector is not large 
enough to encompass two branches of the Fermi surface, 
the attractive potential is broad and shallow enough to 
"put in touch" substantial momentum regions on the dif- 
ferent branches to build up an increased order parameter 
at finite values of cj) away from (j) = 0, 7r/2. 

We also analyzed the case of a CDW instability along 
the (1,1) direction, without finding qualitative differences 
with respect to the (1,0) case. 

Finally we considered an isotropic CDW instability. In 
this latter case the effective attraction is approximately 
given by 



(y/4) 



«2 + |q- qcDVVp 



(8) 



Then we write |q — qcuwP — ((2(2 — cos(ga:) — 
cos((7y))^/^ —qcDwy to guarantee the correct lattice pe- 
riodicity, while keeping the presence of a line of q points 
with large attraction. The result is given in Fig. 5 for the 
case of d-wave symmetry only and for various values of 
qcDW- 
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Fig. 5: Same as in Fig. 4b, but using the potential for an 
isotropic CDW [Eq.(||)]. 

It is apparent that the feature of the local minima at 
(/) = 0, 7r/2 is less pronounced since more q's are effective 
in connecting different points on the Fermi surface. 

To our knowledge there is no report in the literature 
of local minima for the gap at (f> — 0,7r/2. This would 
imply that the qcDW invoked to describe cuprates has 
to be smaller than qp, or the instability should have an 
almost isotropic strucure^ It is worth noting that in the 
Hubbard-Holstein modeld the effective attraction is in- 
deed diverging in an anisotropic way at i5 = 6^^^ ; how- 
ever it has a large value all over the region q ~ qcDW- 
In this case we get an intermediate behavior between the 
anisotropic and isotropic CDW. This case is also very 
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dependent on the interaction couplings and we will not 
show the detailed results aiming to discuss the generic 
features of proximity to CDW. 

In the above discussion of the behavior of the gap near 
CDW we have only partially commented on the s-wave 
solution which also presents rich structures depending on 
ICDW (Figs. 3a and 4a). Indeed we find that the stable 
solution at T = is d-wave even in the proximity to CDW 
similarly to the PS case, provided a sizable U (~ 0.2eV, 
according to our choice) is present. This can be inferred 
from the larger values of the maximum of the d-wave gap 
with respect to the s-wave gap. The explicit computation 
of the condensation energy confirms this expectation in 
all the considered cases. As an example, in Fig. 6, we re- 
port 5F for the set of parameters corresponding to Fig. 4. 
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Fig. 6: s-wave (solid line) and d-wave (dashed line) con- 
densation energy per particle 5F = Fauper — Fnormai as a 
function of the modulus of the CDW wavevector. U — Q.2eV 
and V = 0.3eV. 



For completeness we have also carried out in Appendix 
A a systematic analysis of the relative stability of s- ver- 
sus d-wave solutions in order to find the influence of the 
various parameters (mainly and qcDw) in determin- 
ing the symmetry of the gap function. 



IV. Tc VS DOPING 



The evaluation of a reliable superconducting critical 
temperature is the most difficult issue of a theory deal- 
ing with singular interactions like those generated nearby 
PS and CDW. In the present analysis we start by solving 
Eq.(^ in the linearized form in order to get Tc vs in 
the proximity of PS or CDW (depending on Vg// = Vf^ 



or Vr- 



„ _ YCDW\ 



We report in Fig. 7 the evaluated critical temperature 
for the set of parameters discussed above and for different 
K smaller than the Fermi wave vector kp. 
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Fig. 7: d-wave critical temperatures as a function of the the 
mass parameter close to the PS instability for V = Q.3eV 
(solid line) and to the (anisotropic) CDW instability with 
qcDW = 0.9 and V = 0.3eV (dashed line) and V = 0.45eV 
(dot-dashed line). 

We have verified that the d-wave transition has always 
a substantially larger critical temperature than the s- 
wave (not reported in Fig. 7) both in the proximity of 
PS and CDW. Roughly, the maximum value of the BCS 
gap (at T — 0) gives the order of magnitude of the cor- 
reponding Tc(~ A Max/ '2)- The curves show a strong 
dependence on k^, which at small assumes the form 
Tc oc - logK^. 

The above evaluations of Tc lead to the unphysical re- 
sult that Tc — > oo for S Sc^, Sc^^ , since in these lim- 
its vanishes and the effective potential diverges with 
a non-integrable power. The inclusion of the dynamical 
effects beyond BCS would cutoff the divergency of Tc- 
However, even within the present BCS approach, the sin- 
gularity is cutoff by considering that k^, being the square 
of an inverse correlation length, is vanishing at Sc only 
at T =jJ[)_From the theory of Quantum Critical Points 
(QCP)Lj Ej, it follows that becomes a finite function 
of T for T > const * {S - 5cY^'^'^^^-'^\ where d is the 
dimension, z is the dynamical critical index and j3 is the 
index of the "mass" at T = 0, {5 - 5cY ■ This 

indicates that when PS or CDW occur at finite temper- 
ature, they will occur at 5c{T) < dc{Q) = Sc- Roughly we 
can write 



= Max{a{5 - 5cf ,bT^'^+''-'^^'^) 



(9) 



with a and h model-dependent positive constants, in or- 
der to represent the (much more complex) crossover of 
the actual k'^{5 — 5c, T). Here we will not consider the 
influence of superconductivity on the underlying normal 
phase transition (PS or CDW). This effect would pro- 
vide an additional "mass" contribution expressing the 
stabilizing effect of superconductivity on PS and CDWc3. 
The proper z's are z = 3 for PS and z = 2 for CDW 
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as one can-, see from the inspection of the fluctuation 
propagatord. In d = 2 their values are indeed imma- 
terial since in both cases the above equation reduces 
to = Max{a{S — 6c)^,bT). Dimension two is also 
peculiar because of logarithmic corrections leading to 
k2 ~ a'T(l + 0(logT,log((5 - 6^. This holds in the 
so called classical gaussian regionEa T > {a/b){5 — ScY^ . 

Concerning the index /3, for a QCP one would expect 
fi — 2v — 1, V — 1/2 being thft_classical gaussian index 
of the inverse correlation lengthcil. 

However, PS is a first order transition since symmetry 
does allow the presence of cubic terms (the order param- 
eter being a scalar, density deviation, at g = 0). In the 
parameter space (doping 5 versus bare interaction cou- 
pling g, for instance the electron-phonon coupling in the 
Hubbard-Holstein modelB or thcjifiarest neighbor inter- 
action V in the excitonic modelcjQ) the Maxwell con- 
struction will lead to first order transitions between sta- 
ble phases with > 0. Depending on the models con- 
sidered, the first order transition could end to a criti- 
cal point where PS is second order. Around this point, 
which has a fixed value for both 6 and one would get 

^ {5 — 5cY , i.e. (3 = 2. If the transition is weakly first 
order (i.e. ^ 1) one could also analyse the case (3=1 
for PS. 

The incommensurate CDW is usually of second order 
since cubic terms are not allowed by momentum conser- 
vation. There can be exceptions when qcDW is extremely 
small or when qcDW — 2kF^- For the CDW considered 
in our analysis we assume a second order transition and 
we take (3 = 1. 

From the dependence of = k^{5 — 5c, T) and Fig. 7 
for Tc(«;2) we see that Tc[~ Tc{k^{5 - Sc,T = 0))] 
will rapidly increase by decreasing doping towards dc 
until the value T* ~ {a/bY/'-'^+^~^\5 - 5^)/3^/(d+^-2) 
(T* = {a/b){5 - 5cY in d = 2) is reached at a doping 
5* . Then Tc will slowly approach the value (still of order 
~ T*) which solves the equation 

In the case (3 = 1 and d = 2, this plateau will cover 
all the classical gaussian region Sgl{Tc) < S < S* [GL 
stands for Ginzburg-Landau) . Here 5ql{T) < SdT = 0) 
is the boundary of the tegion where critical fluctua- 
tions start to be relevaniES before arriving to the ac- 
tual critical doping 5c{T) for the normal phase transi- 
tions (PS or CDW), if any. The classical gaussian re- 
gion Sgl{Tc) < S < 6* could be very small and in- 
deed this is the case if the coefficient a is large. In the 
Hubbard-Holstein model we estimate a ~ 10 — 30. On the 
other hand, assuming from Fig. 7 that at optimal doping 
a{S* - Sc) ^ bT^^""" we would get 



'^opt — 0.1 and Kgpj 



(<5* 



JGL) 



{5* 



nlpt/a ~ 0.01. However, be- 



cause of the two dimensional nature of our system, there 
will be a larger region of doping d < 5gl before arriving 
to the normal phase transition. This region will be gov- 
erned by|-q|Uantum fluctuations around the ordered state 
at T = oK 



We suggest that Tc will be a slowly varying function 
of doping in this region, since will mainly depend on 
temperature. Eventually Tc will drop when the normal 
phase transition takes place or the effects of proximity to 
the charge instability are no more effective. In the case 
(3 = 2 \t can also happen that the T = critical point 
is isolated, i.e. no PS instability exists, but for 5 = 5c 
and r = 0. Then we expect a symmetrical behavior at 
5 < 5c with respect to 5 > 5c, and the superconducting 
Tc will fastly decrease as soon as {5 — 5cY > {5* — 5c)'^ 
with 5 < 5c. 

Notice that in discussing Tc versus doping we have as- 
sumed that the main doping dependence is via k? . Indeed 
we have verified that the variations induced by moving 
Ep are less relevant. In Fig. 8 we report the variation of 
Tc by varying doping at fixed and by assuming a dop- 
ing dependence of according to Eq.(^. Specifically, 
Kgpt w 0.1 and = bT for 5 < 5c. Analogous results are 
obtained for AMax versus doping. 
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Fig. 8: Doping dependence of the superconducting criti- 
cal temperature Tc for a fixed value of the mass parameter 
= 0.1. The dotted line is for the case of proximity to PS 
(with V = 0.32eV) and the dashed line is for proximity to a 
CDW instability with \q_cDw\ = 0.9 and V = 0.45eV. In both 
cases the VHS is at 5 = 0.355 and U = 0.2eV. The solid line 
is an estimate of Tc by assuming a doping dependence of 
according to Eq.(^ for d = 2 and (3 = 1 

at fixed K^, the greatest values are obtained for Ep — 
EvHS- Then, a vanishing A.Max and Tc is obtained for 
very high doping {5 > 0.6) or very small (or even neg- 
ative) doping. This analysis is only indicative, since 
the rigid band picture is not valid approaching half- 
filling, where a metal-insulator transition takes place and 
antiferromagnetism appears. Nevertheless it correctly 
describes the rather natural argument that the VHS 
tends to pin the Fermi levelcj, which then requires large 
changes in the filling to be sensibly affected. Then all 
properties which do not directly depend on the rapidly 
varying shape of the Fermi surface (like, e.g. nesting) 
are rather slowly varying functions of the doping. These 
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arguments therefore suggest that strong variations of 
with doping, like those observed in many cuprates, are 
hardly obtained in terms of a dependence on band pa- 
rameters (specifically, tuning the VHS), while they are 
quite natural in the context of proximity to an instabil- 
ity, where doping will control the effective potential itself 
and not only the density of states. This strong variation 
is indeed found when Tc is evaluated by allowing n"^ to be 
doping dependent (solid hne in Fig. 8 to be contrasted to 
the smooth variation of the dotted and dashed curves). 
Notice that our optimal doping value is fixed by proxim- 
ity to the charge instability and not by the proximity to 
the VHS. This agrees with the experimental finding that 
at the maxipaum Tc the VHS is not at the Fermi energy 
but below ito. One can also argue that going away from 
the VHS by decreasing doping could compete with some 
enhancement of k'^{5 — Sc,T) for S < Sc supporting the 
suggestion of a plateau in Tc versus doping for S < dc- 



V. CONCLUSION 

In this paper we have analysed the symmetry of the 
superconducting order parameter in the proximity of PS 
and CDW instabilities. We have found that d-wave pair- 
ing is favored in both cases provided a g-independent re- 
pulsion U > U* is present in the effective interaction 
among quasiparticles. In the CDW case the instabil- 
ity wave vector qcDW should be smaller than the wave 
vectors of the reciprocal lattice, s-wave is realized for 
U < U* depending on qcDW and location of the VHS. 
However, when the s-wave solution is stable, it does not 
present nodes on the Fermi surface both in the proximity 
of PS and CDW in agreement with the finding of Ref.tZl. 
We have considered the band parameters appropriate for 
BiSCCO; however we have found analogous results (not 
reported in the gaper) by using the parameters of YBCO 
as given in RefJj. 

Notice that our static analysis could somewhat over- 
estimate the effects of the repulsive U. However, for the 
range of parameters considered in Sects. II and HI, the 
maximal values U* compatible with a (anisotropic) s- 
wave ground state solution are so small (always smaller 
than O.leV) that it is unplausible that dynamical screen- 
ing would reduce U below these values, thus changing 
our conclusions on the stability of d-wave pairing with 
respect to s-wave pairing. 

We have looked in the Appendix for the factors which 
affects U*, finding that U* increases with increasing 
(i.e. going towards less structured interactions). In the 
case of proximity to CDW instability, an additional pa- 
rameter controlling the d — s-wave interplay is the size 
of the instability wave vector qcDW- Its increasing to 
values of order tt makes the anisotropic density of states 
less relevant for establishing an anisotropic gap. This 
strongly depresses d-wave pairing which indeed disap- 
pears for large qcDW, while it affects less the s-wave 



solution. 

In the context of the model interactions (0) and (0) 
s-wave pairing is "easier" in an electron doped system 
{Ep far away from VHS), specifically if the attraction 
is induced by the proximity to a CDWj-with a sizable 
qcDW- The recent experimental findingl3 on the exis- 
tence of s-wave pairing in overdoped BiSCCO compound 
(if confirmed) is however not easily explained in the con- 
text of proximity to charge instabilities. In particular, 
by overdoping the system approaches the van Hove sin- 
gularity thus enhancing the anisotropy of the gap, which, 
in turn should favor the d-wave symmetry. To reconcile 
this finding with the scenario presented here, one could 
invoke a large reduction of U due to larger screening ef- 
fects, a more three dimensional character of the system, 
or a change of size and direction of qcDW- 

Disorder is an other relevant mechanism which affects 
the relative stability of the d- and s-wave superconduct- 
ing phases: s-wave is more robust in this respect and 
large disorder could favor this latter symmetry with re- 
spect-jto d-wave. This issue was recently considered in 
Ref.CJ and was not of our concern in the present paper. 

Concerning the superconducting critical temperature 
Tc, we have seen that the most important parameter con- 
trolling its value is k^, even though the VHS produce a 
sizable enhancement, which however is effective in a too 
extended region of doping, as seen in Fig. 8. The doping 
dependence of reduces this region. Within the range 
of parameters we have considered (allowing for larger val- 
ues of V within a factor 2), ~ 0.1 ^ 0.2, is enough to 
get Tc ~ lOOK. These values of correspond to cor- 
relation lengths of the order of 2 '--^ 3 lattice constants. 
According to the discussion given in the previous Sec- 
tion, the strong dependence on doping of induces a 
strongly dependent Tc- The maximal Tc is nearby the 
CDW instability for k^{S = 5c, Tc) ~ bTc, with smaller 
values of b favouring higher Tc- Then the maximum Tc 
is not directly related to the proximity to a VHS. 

Various factors affect the coefficient b (like the energy 
scale for the dynamical fluctuations near PS or CDW) 
and its specific estimate depends on the specific model 
leading to charge instabilities. However, we gave argu- 
ments for the existence of a steep increase of Tc with 
doping, just above the QCP for charge instability, fol- 
lowed by plateaus in this dependence and we provided 
rough estimates of this behavior. 

The PS and CDW scenarios have various analo- 
gies withj-Jjlie antiferromagnetic spin fluctuation (AF) 
scenarioBO'EEl. Indeed both approaches assume the prox- 
imity to a QCP. A main claim is that the AF scenario is 
supported by a large number of experimental evidences, 
first of all the actual existence of a AF transition. How- 
ever, notice that the experimental T*^"^ is obtained at 
a doaing far away from 6^^. Moreover, according to 
Ref.O, two remarkable crossover lines can be identified 
in the phase diagram of the superconducting cuprates. 
The first one identifies the doping dependence of the 
maximum in the uniform magnetic susceptibility (Tcr)- 



8 



The second one (T* < Tcr) separates the nuantum crit- 
ical from the quantum disordered regimesu. These two 
curves do not cross, but rathet.converge towards the same 
point [Figs. 2 and 3 in Ref.Elt] nearby the optimal dop- 
ing. It is quite tempting to relate this remarkable point 
with the zero-temperature CDW quantum critical point. 
Indeed, according to our analysis, this point sets up the 
region in doping of maximum attractive interaction and 
should be close to the optimally doped regime. 

Within this scheme, the existence of PS or CDW (once 
long-range forces are taken into account) is not alterna- 
tive to the existence of an AF QCP and the two QCP 
control the behavior of the system at different doping. 
The CDW sets up the maximum critical temperature 
and can constitute the substrate to substain AF fluctu- 
ations far away from the ordered phase, by allowing for 
hole-rich and hole-poor "stripes" . A constructive intee-, 
play between CDW and AF was also suggested in Ref.E2l 
within the analysis of the t-J model. The criticism that 
the experimental J (~ 0.12eM. ~ 0.2t) is lower than the 
value needed for PS (J ~ f)EJ can be overcome by con- 
sidering additional sources for charge. instabilities (for in- 
stance, lecwpling to the lattica3^E3ilj or charge-transfer 
excitongj llil) . With the assumption that the maximal 
Tc is associated to the CDW quantum critical point, the 
quantum disorderd phase will correspond to the region 
between the AF and the CDW quantum critical points. 
Fig. 9 represents this scenario in a schematic way. 
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Fig. 9: Schematic phase diagram for the high temperature 
superconducting cuprates. Both the AF and the CDW quan- 
tum critical point (solid dots) are-phown on the zero tem- 
perature axis. The T« line of Ref.cZI is also indicated by the 
dashed line. The antiferromagnetic (AF), superconducting 
(SC) phases and the spin quantum disordered (SQD) regime 
are also indicated. The crossover line is shaded and phase- 
transition lines are solid. 

According to the above scenario, the occurrence in the 
phase diagram of the high temperature superconducting 
cuprates of both an AF (at low doping) and a CDW (at 
intermediate doping) quantum critical point and the re- 



lated existance of two correlation length, ^ and ^C'^^^ 
opens new possibilities for the interpretation of various 
effects related to the magnetic correlations like, e.g., the 
presence of a spin gap at low and intermediate doping, 
the behavior of the spin-spin correlation length, the dis- 
crepancies epiisting between NMR and neutron scattering 
experimental. In particular, in the spin quantum disor- 
dered regime it is natural to infer the existance of a spin 
gap as a property of spin waves on finite domains of the 
order of 5^^, as also reported in Ref.Ej. 

Within the proposed scenario, owing to the crucial role 
played by CDW fluctuations in determining the physics 
of the superconducting cuprates, it is naturally quite rele- 
vant whether or not an incommensurate CDW symmetry 
breaking occurs. On the other hand the quasi twodimen- 
sional character of the cuprates could play a role in (par- 
tially or totally) suppressing the critical temperature for 
the establishing of a static incommensurate long-range 
CDW order for which a two-component order parameter 
is required. Imlications in this sense are also provided by 
a recent workcj, where the angle-resolved photoemission 
and optical conductivity properties of a model with long- 
range incommensurate CDW order were investigated. A 
better agreement with the observed spectra, was found 
whenever the long-range order was eliminated by the 
(static) superposition of different CDW configurations 
thus mimicking a phase with CDW fluctuations without 
long-range order either due to thermal fluctuations or to 
disorder. This finding suggests that the low dimension- 
ality prevents the establishing in the cuprates of incom- 
mensurate long-range CDW order on a macroscopic scale 
leaving the possibility of static order on a mesoscopic 
scale only or slow dynamical fluctuations. This locally 
(in space and/or time) "ordered" phase should then be 
identified with the magnetic quantum disordered region 
below T, . 

A true CDW order can occur when a suitable match- 
ing between the underlying lattice and the charge fiuctu- 
ations produces commensurability conditions and a con- 
sequent pinning. It was recently proposedtj that such an 
occurrence takes place in La2-x-yNdySrxCu04 at 1/8 
doping. Remarkably, in this systems, the pinned charge 
order and the consequent insulating behaviorj-takes place 
before a long-range AF order is establishedE3. This is 
a clear indication that the freezing of charge fiuctua- 
tions is responsible for the insulating behavior rather 
than the magnetic ordering. Of course this does not 
rule out the possibility of magnetic interactions being 
(co-)responsible for the charge instability, the origin of 
which (magnetic, excitonic, phononic or a combination 
of these mechanisms) is still an open issue. 

The scenario proposed in this paper relates {d-waveX 
superconductivity and anomalous normal state behavioral 
to the proximity to a charge instability. This mecha- 
nism is quite general and is likely expected to work in 
other systems than the high temperature superconduct- 
ing cuprates. A supewionductingr^hase close to both in- 
commensurate SDWES and CDWl^ phases occur in some 
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nearly onedimensional systems, thus suggesting also in 
this case a relation between superconductivity and an 
incommensurate instability. 

Finally we like to mention that superconductivity at 
a fairly large critical temperature (Tc — 31K) occurs 
in Bai-xKxBiOs systems close to a CDW phase. This 
phase at low potassium doping is commensurate and 
responsible for the insulating behavior of the system. 
However, by increasing the doping the system becomes 
metallic and superconducting. Interestingly enough, op- 
tical experimentsEZI show that the feature in the optical 
conductivity related to the CDW gap in the insulating 
phase smoothly and continuously shifts at lower frequen- 
cies upon doping and persist in the mid-infrared region 
in the metallic phase, where it coexists with the Drude 
contribution. It is obviously quite tempting to interpret 
these results as an indication of persistance of CDW in- 
commensurate fluctuations in the metallic phase. If this 
were the case, the scenario proposed here could be of 
relevance for these systems as well. 



APPENDIX A: 

As mentioned in Sects. II and III, a crucial role in 
the relative stability of the d-wave vs anisotropic s-wave 
superconducting phases is played by the local residual re- 
pulsion between the quasiparticles U . In particular, for 
a given parameter set (F, k, qcDw) it is possible to find 
a critical C/* above which s-wave superconductivity be- 
comes unstable with respect to the d-wave phase. This 
quantity is therefore directly related to and provides in- 
formation on the relative robustness of the two supercon- 
ducting phases: The larger [/* is and the more difhcult 
it is to spoil s-wave superconductivity. In order to filter 
this information from the absolute strength of supercon- 
ductivity for a given parameter set, we normalize the [/* 
with {Vattrigx^Qy)), the average value (on the Brillouin 
zone) of the attractive part of the effective potentials in 
Eqs. d) and (0). 

We carried out an analysis of the normalized critical lo- 
cal repulsion U* / (Vattriqx, %)) as a function of the mass 
parameter for the PS instability, and as a function of 
the wavelength for the CDW instabilities with wavevec- 
tors in the (1,0) and (1,1) directions and for the isotropic 
CDW. Table I displays the results of this analysis. Two 
clear tendencies are found. First of all the s-wave solu- 
tion is made more stable {U* / {Vattr{qx , qy)) increases) by 
increasing the k^. This is because in this way the effec- 
tive potential becomes less sharply structured in momen- 
tum space, the attractive well being more shallow. As a 
consequence, the superconducting gap anisotropy is less 
pronounced (there is a weaker momentum decoupling) 
making in turn comparatively less favorable to create a 
line of nodes in the order parameter to produce d-wave 
superconductivity. 

The second clear tendency is the increase of 



U* /{Vattr{qx,qy)) by increasing the momentum of the 
CDW instability. Increasing |qcDw| favors the couphng 
of momenta, which are rather distant from each other, 
thus reducing the effect of the anisotropic density of 
states. Consequently, the gap anisotropy is smeared and 
also in this case the lines of nodes of the c?-wave solution 
are made comparatively less favorable. Notice also that 
this latter effect is more pronounced for the instability in 
the (1,1) direction because the phase of the order param- 
eter in the d-wave superconductivity is hardly compati- 
ble with an attractive potential strongly scattering states 
close to the (tt, 0)-(0, tt) points of the Brillouin zonecj. 

We finally like to point out that fixing {Vattriqx, qy)) in 
order to obtain (for the pure d-wave) a critical tempera- 
ture around lOOK, in all cases U* turns out to be quite 
small (a few hundreds of eV). 



10 



PS 




CDW (1,0) 




CDW (1,1) 




CDWiso 




K 




qcDw 


U*/{V) 


qcDw 


U*/{V) 


qcDw 


U*/{V) 


7r/16 


0.35 


0.9 


0.40 


0.9 


0.42 


0.9 


0.29 


7r/8 


0.42 


1.2 


0.45 


1.2 


0.49 


1.2 


0.42 


7r/4 


0.51 


1.5 


0.52 


1.5 


0.59 


1.5 


0.59 






1.8 


0.56 


1.8 


0.74 







TABLE I. Dependence of the normalized critical U* on 
the mass parameter n close to a PS and on the modulus of 
the instability wavevector close to the CDW instability (for 
= 7r/16). 
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